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Abstract 

Modified Rayleigh conjecture (MRC) in scattering theory is proposed and justi- 
fied. MRC allows one to develop numerical algorithms for solving direct scattering 
problems related to acoustic wave scattering by soft and hard obstacles of arbitrary 
shapes. It gives an error estimate for solving the direct scattering problem. It sug- 
gests a mimcrical method for finding the shape of a starshaped obstacle from the 
scattering data. 



1 Introduction 

Consider a bounded domain D C M", n = 3 with a boundary S. The exterior domain 
is D' — W?\D. Assume that S is smooth and starshaped, that is, its equation can be 
written as 

r = /(«), (1.1) 

where a G 5^ is a unit vector and S"^ denotes the unit sphere in E?. Smoothness of S is 
used in (4.6) below. For solving the direct scattering problem by the method described 
in the beginning of Section 2, the boundary S can be Lipschitz. The acoustic wave 
scattering problem by a soft obstacle D consists in finding the (unique) solution to the 
problem (1.2)-(1.3): 

(V^ + A;2) = in D', = on S, (1.2) 

u — uq -\- Aia\ a) h o - ) , r := b| ^ oo, «':=—. (1.3) 

r \r J r 

*key words: Rayleigh hypothesis, scattering by obstacles, inverse scattering 
■I^AMS subject classification: 35R30 



Here Uq := e*'^"'^ is the incident field, A{a', a) is called the scattering amplitude, its 
k-dependence is not shown, > is the wavenumber. Denote 

Ai{a):= / A{a',a)Yi{a')da', (1.4) 

where Y£{a) are the orthonormal spherical harmonics, = Yim, —i<m<i. Let he{r) 
be the spherical Hankel functions, normalized so that hi^r) ~ as r +00. Let the 
ball Bn := {x : |x| < R} contain D. 

In the region r > R the solution to (1.2) - (1.3) is: 

00 

m(x,«) =e''="-" + 5^A,(«)^,, iJe:=Y,{a')h,{kr), r > R, a' = ^, (1.5) 
e=o 

summation includes summation with respect to m, —£ < m < £, and Ae{a) are defined 
in (1.4). 

Rayleigh conjecture (RC): the series (1.5) converges up to the boundary S (originally 
RC dealt with periodic structures, gratings). This conjecture is wrong P, [3], [4]. For 
example, if n = 2 and D is an ellipse, then the series analogous to (1.5) converges in the 
region r > a, where 2a is the distance between the foci of the ellipse j]P . In the engineering 
literature there are numerical algorithms, based on the Rayleigh conjecture. Our aim is 
to give a formulation of a modified Rayleigh conjecture (MRC) which is correct and can 
be used in numerical solution of the direct and inverse scattering problems. We discuss 
the Dirichlet condition but similar argument is applicable to the Neumann boundary 
condition, corresponding to acoustically hard obstacles. 

Fix e > 0, an arbitrary small number. 

Lemma 1.1. There exist L = L{e) and q = Q(e) such that 

L{e) 

||mo + X^Q(e)V'«IU2(5) < e- (1-6) 
£=0 

If (1.6) and the boundary condition (1.2) hold, then 

L{e) 

\\ve - v\\l2(s) < e, V, ■.= '^ce{e)i)i. (1.7) 

1=0 

Lemma 1.2. If (1.7) holds then 

\\v^-v\\ = 0(e) e^O, (1.8) 

where || ■ || := || ■ \ \hj^^{D') + II " I |l2(d';(i+|x|)-^)) ■J > I, m> is an arbitrary integer, and 
is the Sobolev space. 
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In particular, (1.8) implies 

IK-^||L2(5K) = 0(e) e^O. (1.9) 

Lemma 1.3. One has: 

ce{e) ^ Ae{a)Wi, e ^ 0. (1.10) 

The modified Rayleigh conjecture (MRC) is formulated as a theorem, which follows 
from the above three lemmas: 

Theorem 1 (MRC): For an arbitrary small e > there exist L(e) and Q(e),0 < 
e < L{e), such that (1.6), (1.8) and (1.10) hold. 

The difference between RC and MRC is: (1.7) does not hold if one replaces by 
S^=o ^i{^)'^tj let L ^ oo (instead of letting e —>■ 0). 

For the Neumann boundary condition one minimizes 1 1 9b^o+J2e=o ^ei>e] 1 1^2(5) with re- 
spect to cg. Analogs of Lemmas 1.1-1.3 are vahd and their proofs are essentially the 
same. 

In Section 2 we discuss the usage of MRC in solving the direct scattering problem, in 
Section 3 its usage in solving the inverse scattering problem, and in Section 4 proofs are 
given. 



2 Direct scattering problem and MRC. 

The direct problem consists in finding the scattered field v given S and uq. To solve 
it using MRC, fix a small e > and find L(e) and Q(e) such that (1.6) holds. This is 
possible by Lemma 1.1 and can be done numerically by minimizing | |mo+X1o ^ii^^WL^is) '■ = 

4>{ci, ,cl). If the minimum of is larger than e, then increase L and repeat the 

minimization. Lemma 1.1 guarantees the existence of such L and q that the minimum is 
less than e. Choose the smallest L for which this happens and define :— X]t=o ^ei^ei^)- 
Then is the approximate solution to the direct scattering problem with the accuracy 
0(e) in the norm || ■ || by Lemma 1.2. 

In [6] representations of v and are proposed, which greatly simplified minimization 
of 0. Namely, let solve problem 

(V^ + k^) = in D', = on S, (2.1) 

and satisfies the radiation condition. Here {fe}e>o is an arbitrary orthonormal basis 
oiL\S). Denote 

00 

v{x) := Ci'^eix), u{x) := uo + v{x), q := {-uq, fi)L^{s)- (2.2) 

The series (2.2) on jS" is a Fourier series which converges in LF'{S). It converges 
pointwise in D' by the argument given in the proof of Lemma 1.2. A possible choice of 
fi for star-shaped S is fi = Yij ^Jw where w := dS/da. Here dS and da are respectively 
the elements of the surface areas of the surface S and of the unit sphere S^. 
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3 Inverse scattering problem and MRC. 



Inverse obstacle scattering problems lOSPa) and lOSPb) consist of finding S and the 
boundary condition on S from the knowledge of: 

lOSPa): the scattering data A{a', a, ko) for all a', a e S^, k — ko > being fixed, 

or, 

lOSPb): A{a', a^, k), known for all a' G S'^ and all A: > 0, a = cio G S'^ being fixed. 

Uniqueness of the solution to lOSPa) is proved by A.G.Ramm (1985) for the Dirichlet, 
Neumann and Robin boundary conditions, and of lOSPb) by M.Schiffer (1964), who as- 
sumed a priori the Dirichlet boundary condition. The proofs are given in [4]. A.G.Ramm 
has also proved that not only S but the boundary condition as well is uniquely defined 
by the above data in both cases, and gave stability estimates for the solution to lOSP [9]. 
Later he gave a different method of proof of the uniqueness theorems for these problems 
which covered the rough boundaries (Lipschitz and much rougher boundaries: the ones 
with finite perimeter [8], see also [10]. In [11] the uniqueness theorem for the solution 
of inverse scattering problem is proved for a wide class of transmission problems. It is 
proved that not only the discontinuity surfaces of the refraction coefficient but also the 
cocflicicnt itself inside the body and the boundary conditions across these surfaces are 
uniquely determined by the fixed-frequency scattering data. For any strictly convex, 
smooth, reflecting obstacle D analytical formulas for finding S from the high-frequency 
asymptotics of the scattering amplitude are proposed by A.G.Ramm, who gave error 
estimates of his inversion formula also [4] . The uniqueness theorems in the above refer- 
ences hold if the scattering data are given not for all a', a e S'^, but only for a' and a 
in arbitrary small sohd angles, i.e., in arbitrary small open subsets of S'^. The inverse 
scattering problem with the data a' & S^, k — ko and a — ao being fixed, is open. If a 
priori one knows that D is sufficiently small, so that /cq > is not a Dirichlet eigenvalue 
of the laplacian in D, then uniqueness of the solution with the above non-overdetermined 
data holds (by the usual argument [4]). There are many parameter-fitting schemes for 
solving lOSP, [13], see also [5]. 

Let us describe a new such scheme, based on MRC, its idea is similar to the one 
in [7]. Suppose that the scattered field v is observed on a sphere Sr. Calculate q 
{v,Yi)i2(^g2-j/hi{kR). If V is known exactly, then q = Ag{a). If vg are noisy data, ||f — 
vs\\l^(Sr) ^ ^) then q = C£s- Choose some L, say L = 5, and find r = r(a') as a positive 
root of the equation uq + vl := e**^"'"''' -|- '^^^QCis'4'i{kr,a') := p{r,a',a,k) = 0. Here 
a' and k > are fixed, and we are looking for the root r = r{a') which is positive and 
stable under changes of k and a. In practice equation p{r, a',a,k) — may have no such 
root, the root may have small imaginary part. If for the chosen L such a root (that is, a 
root which is positive, or has a small imaginary part, and stable with respect to changes 
of k and a) is not found, then increase L, and/or decrease L, and repeat the search of 
the root. Stop the search at a smallest L for which such a root is found. The MRC 
justifies this method: for a suitable L the function p(r, a', a, k) is approximately equals 
zero on S, that is, for r = r{a'), and this r{a') does not depend on k and a. Moreover, 
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by the uniqueness theorem for lOSPa) and lOSPb) there is only one such r = r(a'). 
Numerically one expects to find a root of the equation p{r, a', k) = which is close to 
positive semiaxis r > and stable with respect to changes of k and a. 

If one uses the above scheme for solving the inverse scattering problem for an acous- 
tically hard body (the Neumann boundary condition on S), then one gets not a tran- 
scendental equation p{r,a',a,k) = for finding the equation of S', r = r{a'), but a 
differential equation for r = r(a'), which comes from the equation ^^^''^'"'^^ = at 
r = r(a'). One has to write the normal derivative on 5* in spherical coordinates and then 
substitute r = r(a') into the result to get a differential equation for the unknown function 
r = r{a'). For example, if n = 2 (the two-dimensional case), then the role of a' plays the 
polar angle if' and the equation for r = r{ip') takes the form ^ = {r'^^/ -^)\j.=r(ip')- 

4 Proofs. 

Proof of Lemma 1.1. This Lemma follows from the results in [4], (p. 162, Lemma 1). 
Proof of Lemma 1.2. By Green's formula one has 



Ve{x) = - / Ve{s)GN{x,s)ds, \\ve{s) +mo||l2(5) < e, (4.1) 
Js 

where G is the Dirichlet Green's function of the Laplacian in D': 

(V^ + k^)G = -6{x - y) in D\ G = on S, (4.2) 
dG 



lim / —-- ikG ds = (4.3) 

r^<^Jl^l=r d\x\ 

From (4.1) one gets (1-3) with i7/^^(D')— norm immediately by the Cauchy inequality, 
and with the weighted norm from the estimate 

\Gn{x,s)\<—^, \x\>R, (4.4) 
1 + \x\ 

and from local elliptic estimates for '■= — v, which imply that 

\\we\\L^{Bii\D) < ce. (4.5) 

Let us recall the elliptic estimate we use. Let D'j^ := Br\D, Sr be the boundary of 
Bn, and choose R such that /c^ is not a Dirichlet eigenvalue of —A in D'j^. The elliptic 
estimate we have used is ([2], p. 189): 

\\We\\H"^{D'^) < C[||(A + k'^)We\\H^-^D'^) + | 1 1 //™-0.5 (5^) + 1 1 U^e | U'"-0-5(5)] • (4.6) 

Take m = 0.5 in (4.7), use the equation {A + k'^)w^ = in D', the estimate | |we| |_H-"i(s'fl) = 
0(e), proved above, the estimate | |_h-o(5) = 0(e), and get (1.8). Lemma 1.2 is proved. 
□ 
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Proof of Lemma 1.3. Lemma 1.2 yields convergence of to v in the norm || ■ ||. 
In particular, \\v^ — v\\l2(^Sr) — as e — 0. On Sji one has v = Y^'^^ Ai{a)ipi and 
= ^ftoCeipe- Multiply Ve{R,Oi') — v{R,a') by Yi{a'), integrate over 5"^ and then let 
e ^ 0. The result is (1.10). □ 
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